Bose-Einstein condensation transition studies for atoms confined in 
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Multiply-connected traps for cold, neutral atoms fix vortex cores of quantum gases. Laguerre- 
Gaussian laser modes are ideal for such traps due to their phase stability. We report theoretical 
calculations of the Bose-Einstein condensation transition properties and thermal characteristics of 
neutral atoms trapped in multiply connected geometries formed by Laguerre-Gaussian (LGp) beams. 
Specifically, we consider atoms confined to the anti-node of a LGj laser mode detuned to the red 
of an atomic resonance frequency, and those confined in the node of a blue-detuned LGj beam. 
We compare the results of using the full potential to those approximating the potential minimum 
with a simple harmonic oscillator potential. We find that deviations between calculations of the full 
potential and the simple harmonic oscillator can be up to 3% — 8% for trap parameters consistent 
with typical experiments. 
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Bose-Einstein condensation (BEG) has led to a wealth 
of new physics, from applications in inertial measure- 
ments to fundamental studies of statistical mechanical 
phenomena and superfluidity jli]. Among the first experi- 
ments probing BEG were those that explored the thermal 
properties and transition characteristics, with particular 
emphasis on the effect of trap geometry 0, y] . An im- 
portant avenue of investigation in the connection between 
condensed matter (super-fluidity) and BEG [J| of atomic 
gases is the experimental [^-Q and the theoretical [1, Q 
studies of vortices. It was recognized early that a trap 
potential with the ability to pin a vortex core would be an 
advance in the study and application of these quantized 
rotations in BEG [l3l • Toward this end, researchers have 
developed novel multiply-connected trap geometries. A 
toroidal potential locally fixes a vortex state, thus en- 
couraging rotational stability for more precise measure- 
ments and may help facilitate the development of devices 
utilizing vortices, such as the atomic squid detector [lli| 
and the ultra-stable gyroscope pJi]. Recent progress in 
toroidal potentials has been seen in the confinement of 
cold atoms [Tj], BEG [Tj, [l^, and creation of vortex 
states p^ . 

An important method for creating toroidal geome- 
tries uses higher-order Laguerre-Gaussian laser modes. 
The Laguerre-Gaussian beams, as a set of solutions to 
Maxwell's equations T6], represent stable modes of laser 
propagation. The radial electric field is proportional to 
an associated Laguerre polynomial, L' and a Gaussian 
function. The electric field of LG' laser modes has an az- 
imuthal winding phase given by e"*'*^ [la|, and the p+l 
radial intensity nodes provide a variety of multiply con- 
nected geometries for vortex studies. LG' photons have 
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a quantized orbital angular momentum (QAM) of Ih per 
photon. This property can also be exploited to create 
vortices throiigh the coupling of photon-matter QAM, as 
proposed by [13, [l3 and demonstrated by [ifl [3 ■ 

We explore the thermal properties (the transition tem- 
perature, the population of atoms in the ground state, 
and the specific heat) of a Bose gas trapped in a LG^ 
dipole trap, utilizing the complete LG' potential. These 
results are compared to those where the confining poten- 
tial minimum has been approximated as a simple har- 
monic oscillator. Previously, the simple harmonic oscil- 
lator approximation has been used in calculations to de- 
termine the ground state energies, density profiles, and 
the transition temperature in LG beams [20l422l |. 

Gonsider an atom in an inhomogeneous laser field 
whose angular frequency, lo, is near enough to a reso- 
nant angular frequency, wq, that the coupling between 
any other pair of states can be ignored. The laser detun- 
ing, A = a; — Wo, is sufficiently large that the probability 
of photon absorption is negligible. The atom experiences 
an attractive force toward regions of high laser intensity 
when the detuning is negative (A < 0) and a repulsive 
force away from regions of high laser intensity when the 
detuning is positive (A > 0). 

We investigate two physically realizable toroidal trap- 
ping geometries formed by the lowest two orders in a LG' 
beam (with non-zero angular momentum): the LGj red- 
detuned trap and the LG} blue-detuned trap. The in- 
tensity profile for the LG J mode is given in Figure [1] (a) . 
For A < 0, this shape satisfies the desired toroidal geom- 
etry with atoms confined to regions of high laser inten- 
sity. Blue detuned dipole traps are preferred for applica- 
tions where external perturbation needs to be minimized. 
Atoms spend the majority of their time in regions of low- 
est intensity. The LG} beam satisfies this condition, and 
a sample profile is shown in Figure [1] (b) . 
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FIG. 1. (Color online) The transverse intensity profile for (a) 
the LGo mode and (b) the LGi mode. 



The potential ener gy o f an atom in the presence of a 



LGp mode is given by [21'], 
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where Pq is the total laser power, T = 27r x 6.1 MHz is 
the natural line width for ^ Rb, and w is the beam waist. 
The resonant saturation intensity is Isat = tt/ic/SA^t, 
where A is the resonance wavelength and r is the lifetime 
of the resonant state. We work in cylindrical coordinates 
where the dipole potential only confines in the transverse 
r-dimension. To confine along the direction of laser prop- 
agation (z-dimension), we will assume a harmonic poten- 
tial. One dimensional confinement can be constructed via 
an anisotropic magnetic trap whose trapping potential in 
the radial dimension is small enough to be ignored com- 
pared to the laser potential [23j, ]2M ■ With the inclusion 
of the 1-D harmonic confinement together with the radial 
dipole potential in Equation ([T]), the confining potential 
for both laser modes has the form, 
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where ujz is the trapping frequency in the z-direction and 
TO is the mass of the confined atom. The constant V is 



FIG. 2. (Color online) The radial potential energy (in units 
of euo) of an atom in specific LGp laser modes as a func- 
tion of distance from the center of the laser beam (in units 
of aho)- In both plots, the solid line corresponds to the full 
LGp dipole potential, and the dashed line represents the sec- 
ond order expansion about the minima of the potential (SHO 
approximation) . 



the product of three physical quantities, 
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The first factor is an energy scale set by the energy width 
of the excited state. The second factor is the ratio of an 
intensity term (/q = Pq/ttw'^) to the saturation intensity. 
The last factor is a ratio of the line width to the laser 
detuning. These factors combine to set the scale of the 
trap strength. 

The 1-D radial cross-section of both potentials in 
Equation ^ are plotted in Figure[2] The calculations are 
done in di mensionl ess units. All the lengths are scaled 
by Qho — ^/fi/mujz, the harmonic oscillator length of the 
simple harmonic confining potential in the z-direction. 
All the energies are scaled by Cho = f>^z/2, the en- 
ergy of the corresponding 1-D ground state in the z- 
direction. Dimensionless quantities are identified with a 
tilde {'^) set above the term: f = r/aho, z = z/aho, w — 
w/ufio, e = e/eho, and v = V/eho. Making these substi- 
tutions into the potentials. Equation ([2]) becomes. 
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We explore a few of the thermal properties of a cold 
Bose gas confined in a toroidal potential. One such prop- 
erty is the BEG transition temperature, Tc- We calcu- 
late Tc following the procedure outlined in [23 ■ In the 
limit when N ^ 00 and ksT is much larger than the 
energy spacing of the confining potential, then 



N = N, 
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where Nq is the number of atoms in the ground 
state, p(e) is the density of states, and n(e) is the usual 
Bose distribution occupation number l/{e^^'^~f^'> — 1), 
where /3 = l/fc^T and fi is the chemical potential. The 
density of states is 
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where V*{e) is the classical spatial volume spanned by a 
particle with energy e. 

As the phase space of the system decreases, the ground 
state population remains unoccupied until the chemical 
potential approaches zero. Thus, the transition temper- 
ature, Tc, can be found by setting fj. — and Ng = 0. 
The heat capacity, C'(T), is piecewise continuous, but a 
discontinuity at Tc is a signature of the BEG phase tran- 
sition. The heat capacity is C{T) = dE{T)/dT where 
the total energy of the system is 
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Equations dH) through ([6]) can be solved analytically 
when the potential is expanded up to the second order 
about the potential minimum. This is the simple har- 
monic oscillation (SHO) approximation used in |20l - |22| . 
For our system and notation, the confining potentials are 
given by 
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in the SHO approximation. Defining {pp)ho 
specific LGp potential, the densities of states are 
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From here, we carry out the integrations in Equation (|4]) 
to find the BEG transition temperature, {Tp)ho, 
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where the term C(5/2) is the Riemann Zeta function. 
Galculation of Tc using the full potential given 



by Equation ([3]) must be done numerically. For a 
given w and v, the density of states (Equation ([5])) is nu- 
merically integrated for a discrete set of energies, e^. The 
numerical integration is done in Mathematica [26| and 
wc provide our source code online [27|. Twenty values 
of the density of states are calculated covering an en- 
ergy range ei = to emax- The value of emax is de- 
termined such that most of the Bose-Einstein distribu- 
tion is accounted for, while being certain that the energy 
range sufficiently maps out the most probable portions 
of the trap occupied by the particle. To do this, we find 
the location of the critical temperature of the SHO ap- 
proximation with respect to the potential barrier of the 
corresponding LG^ mode. The energy range is adjusted 
to span this region. The limits of integration for Equa- 
tion ^ are determined by solving Equation (jS]) for the 
classical volume accessible by a particle of energy e^ . For 
the radial classical turning points, we use the Newton- 
Raphson numerical technique. The resulting set of points 
of the density of states, Pi(ei), are fit to the analytical ex- 
pression. 
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where the factors cr and rj are the fitting parameters. 
This model has a general form that is characteristic for 
many systems. For example, a Bose gas with no exter- 
nal potential results in 77 = 1/2, one confined by a 3-D 
isotropic oscillator gives rj — 2, and both harmonic ap- 
proximations to the LG beams given in Equations ^ and 
([5]) result in ?7 = 3/2. We substitute the fitted function 
into Equation Q and analytically integrate to find Tc for 
a given N (in our case, we choose a value of A^ = 10®). 
Finally, we analytically integrate Equation ^ to find the 
total energy, E{T), and then the heat capacity. 

Figures [3] through [7] show the expected BEG tran- 
sition temperature, ground state fraction, and heat ca- 
pacity for a ®^Rb gas confined in the LG' modes. We 
also show the corresponding harmonic approximations 
and compare with the exact calculation for different laser 
beam parameters. Figure [3] shows Tc in the LGj (up- 
per line) and the LGj (lower line) modes as a function 
of the trap depth while maintaining equal trap frequen- 
cies in the axial and radial directions. The trap depth 
is the height of the lowest confining barrier of the ra- 
dial potential (Figure [5]). The axial trap frequency is 
set to ujzjI-K = 100 Hz, typical of magnetic trap ex- 
periments [1]. The trap depth is changed by adjust- 
ing V, equivalent to adjusting the power or detuning of 
the trapping laser. At the same time, w is adjusted 
so the radial trap frequency is maintained at ujr = i^z- 
The critical temperature is calculated using the full po- 
tential (solid line) and using the SHO approximation of 
the potential well (dotted line). Note that Tc decreases 
with the trap depth in this figure because of the con- 
straint ujr = Wz- As the laser power increases, the beam 
waist must increase to maintain constant ujr- This in- 
creases the trap volume, and decreases Tc- Figure [3] (b) 
shows the fractional error of the SHO approximations 
compared to the exact calculation. The uncertainty is 
largest when the transition temperature approaches the 
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FIG. 3. (Color online) (a) Tc calculations using both the full 
potential (solid curve) and the SHO approximation (dashed 
curve) as a function of trap depth, where the radial trap fre- 
quency is held equal to the axial trap frequency, uiz/'ii^ ~ 
lOOHz. (b) The factional error of the SHO calculation. 



trap depth. The classical volume accessed by particles 
includes a larger fraction of the potential well where the 
harmonic approximation is no longer valid. We calculate 
fractional differences from 7.6% to 0.023% over a trap 
depth of 0.1 fiK to 10 fjK. 

In contrast to Figure |31 Figure 0] considers a sym- 
metric case where the beam waist is held constant while 
increasing the trap depth. This corresponds to the trap 
volume decreasing as the trap depth increases. The 
critical temperature in the LGq (upper line) mode and 
the LGj (lower line) mode is plotted as a function of the 
laser power. We fix the beam waist at w = 42 //m for the 
LGj mode, and w = 58 /j,m for the LGj mode. The laser 
detuning is set to |A|/27r = 1000 GHz for both modes. 
The laser power is varied over a range of 0.52 mW— 1.3 W 
for the LGj mode, and 1.3 mW - 5.9 W. This has the ef- 
fect of increasing the trap depth and Ur ■ To keep the trap 
symmetric, we allow Uz to change such that ujz = uir- An 
increase in the trap depth has an effect of increasing Tc- 
For this case, we also make the SHO approximation to 
the confining potential. The fractional differences are 
from 4.5% to 0.017%. Like in Figure |31 the largest error 
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FIG. 4. (Color online) (a) Tc as a function of the laser 
power. The beam parameters are Wr = 50 /im and A/27r = 
1000 GHz. The comparison is made using both the full poten- 
tial (solid curve) and the SHO approximation (dashed curve) 
while keeping the trapping beam waist fixed and simultane- 
ously increasing the laser power. In addition, the axial har- 
monic trapping frequency is adjusted such that the trap re- 
mains symmetric at all laser powers, (b) The fractional error 
of the SHO calculation. 



occurs when the trap depth and the transition tempera- 
ture are comparable. As a larger fraction of the particles 
populate larger regions of the potential, the particles oc- 
cupy regions where the potential becomes less harmonic. 
We replace the constant, a in Equation (fTO|) . with the 
polynomial a{e) = ctq + trie -I- a'2e^ • • • , where the min- 
imum number of terms were kept such that the result 
converged. 

Figure [5] considers the effect on Tc when the trap- 
ping frequencies in the two trapping dimensions (r and 
z) are allowed to vary with respect to each other. The 
symmetric case (aJr/27r = 100 Hz on the horizontal axis 
in Figure[5]) corresponds to a laser power of Pq = 10 mW, 
a beam waist w = 50 /im, and a detuning of |A|/27r = 
1000 GHz. These values are concurrent with typical har- 
monic oscillator lengths of aho = 1 IJ-^ca reported by [l| . 
We vary the radial trap frequency by fixing the waist of 
the trapping laser and increasing or decreasing the laser 



300 



250 
c 200 



k" 



150 



100 




100 150 200 300 
Radial Trap Frequency (Hz) 




100 150 200 300 
Radial Trap Frequency (Hz) 

FIG. 5. (Color online) (a) Tc calculations using both the full 
potential (solid curve) and the SHO approximation (dashed 
curve) while keeping the trapping beam waist fixed and si- 
multaneously increasing the laser power. Thus, the effects of 
adjusting the radial trapping frequency with respect to the 
axial frequency. Again, the axial frequency is fixed at lOOHz. 
(b) The factional error of the SHO calculation. 



power. In terms of the quantities defined in this paper, 
this is the same as keeping w fixed while increasing v. Un- 
like the results in Figure [3l here the confinement in the z- 
dimension is held constant leading to asymmetric traps. 
We calculate T^ assuming the SHO approximations to 
the confining potential. For radial frequencies uj^ < lOz, 
the atoms experience weaker confinement in the radial 
dimension, and the trapped atoms access regions of the 
trap where the potential is less harmonic. The fractional 
error for the LGj (top line) and the LG} (bottom line) 
are plotted in Figure [5] (b). For a range of frequen- 
cies spanning 40 — 400 Hz, we find fractional differences 
of 2.7-0.07%. 

We have also performed similar calculations for atoms 
trapped in a red detuned LG} laser beam. This confining 
potential is more complicated since it provides concentric 
multiply-connected traps. These concentric traps have 
different trap depths and asymmetries. However, the re- 
sults are consistent with the calculations of the blue de- 
tuned LG} and red detuned LGj modes - BEG transition 
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FIG. 6. (Color online) (a) Plot of the fraction of particles in 
the ground state as a function of the scaled temperature. Plot- 
ted are the SHO approximation to the LG potentials (black 
dashed), the LGi mode (blue solid), the LGq mode (red dot- 
dashed), and the 3-D isotropic harmonic oscillator (green dot- 
ted), (b) The absolute difference in LGp fractions and the 
fraction due to the SHO approximation 



characteristics calculated from the full potential deviate 
from those using a SHO approximation when the clas- 
sical volume of the potential accessed by particles when 
the system is near Tc includes regions where the harmonic 
approximation to the potential breaks down. 

From Equation ^ the fraction of atoms in the ground 
state is 
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where /3c — ^/kBTc- Figure[6](a) shows Nq/N as a func- 
tion of the scaled temperature, T/Tc for the LGq mode 
dipole trap (dot-dashed line), the LG} mode dipole trap 
(solid line), the harmonic approximation to a LG' mode 
(dashed line), and the 3-D SHO (dotted line). The trap 
parameters for the LGj mode consist of a laser power of 
Pq = 0.16 mW and a beam waist oi w — 31 /im. The 
trap parameters for the LG} mode are Pq — 1.3 mW 
and a beam waist of w = 58 mW. Both modes have a 
laser detuning of \A\/2tt = 1000 GHz and have the same 
trapping frequency, w^ = a;,, = 27r x 100 Hz. Note that 
given the form of the density of states (Equations (j8]) 




FIG. 7. (Color online) Plot of the heat capacity per particle 
in units of kt for the LGq (red) , the LGi (blue) , and the har- 
monic approximation (black) . The heat capacity is a function 
of the temperature, where T is in units of Tc- The thermody- 
namic discontinuity at T — Tc indicates a phase transition. 



and (ITU|) ). Equation (|lip does not depend on the coeffi- 
cient, a, but only on the exponent, rj. Therefore, Nq/N is 
exactly the same for the SHO approximation to either po- 
tential. 

The 3-D SHO is included in Figure H] (a) to give a qual- 
itative perspective. The errors associated with making 
the harmonic approximation of the LG modes are sim- 
ilar to the differences between the multiply connected 
toroidal geometries and the singly-connected 3-D SHO 
potential. Figure|6](b) reports the absolute difference be- 
tween the two LG modes and their respective harmonic 
approximations. The harmonic approximation error can 
reach as high as 3.7% for the LGj and 1.8% for the LG}. 
Such errors are on the same order as the effects associ- 
ated with finite system sizes {N — > cx) approximation) 
experimentally shown in reference [3|. The trap depth 
is 0.5 yuK, the beam waist for the LGj is 31 fim, the 
beams waist for the LG} is 58 /im, and the laser detun- 



ing is |A|/27r = 1000 GHz. A precise knowledge of the 
fraction of trapped atoms in the ground state may be 
necessary for many of the applications of LG laser modes 
as dipole traps. 

Figure [7] shows the heat capacity for samples trapped 
in LGj (upper curve) and LG} (middle curve) with a re- 
spective harmonic approximation (lower dashed curve). 
Like Figure [Bl the harmonic approximation can be rep- 
resented by a single curve because the scaling of the heat 
capacity and temperature allow for a dependence on a 
single parameter, rj in Equation (jlOp . The laser param- 
eters correspond to situations where the confined atoms 
occupy regions of the trap that deviate from a harmonic 
description. The disagreement between the LG poten- 
tials and the harmonic approximations can be as large 
as 24% around Tc. 

Toroidal traps are ideal geometries for employment in 
BEC gyroscope and vortex experiments. Matter-wave in- 
terferometers may offer better sensitivity than traditional 
optical interferometers. Vortices provide a key testing 
ground for superfluid behavior studies. For these precise 
applications, BEC characteristics should be well under- 
stood. In particular, knowledge of the fraction of trapped 
atoms in the condensate is important for these applica- 
tions. In addition, the wave functions of condensates 
confined in LG traps may need to take into account the 
full trap geometry. 

We calculate the thermal properties of a Bose gas con- 
fined by a LGp laser dipole trap: the critical transition 
temperature for a wide range of experimental parameters, 
the condensate fraction, and the heat capacity. We also 
compute the thermal properties using a SHO approxima- 
tion to the potential minimum and compare to the exact 
results. Depending on the precision required, there ex- 
ists a regime in which the thermal properties of atoms 
confined to LG^ dipole traps contain non-negligible de- 
viations. When the critical temperature is on the order 
of the trap depth (a depth of 0.1 /iK), it is predicted to 
be too high by errors as large as 8%. In this regime, 
the number of atoms in the ground state is underesti- 
mated for temperatures below the critical temperature. 
We also find that the heat capacity contains large devi- 
ations around the BEC transition temperature. These 
corrections are on the order of well known effects (such 
as finite size corrections). 

This work is supported by the Research Corporation, 
Digital Optics Corporation, and The University of Okla- 
homa. 
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